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ABSTRACT 

The first detailed VLA mapping of the neutral hydrogen distribution in the isolated, 
edge-on spiral galaxy NGC 3044 is presented. Physical parameters such as M HI , M T , 
etc. determined for this galaxy are typical for galaxies of its morphological class (SBc). 
We have modelled the HI spectra in order to derive its global density and velocity 
distributions. An HI scale height of 420 h~ l pc is thus found. This can be compared 
to the impressive radio continuum halo, previously found to extend to 8 kpc above the 
midplane. 

The present study reveals an asymmetry in the HI distribution as well as numerous 
high-latitude HI structures at various galactocentric radii. The approaching (northwest) 
side of the galaxy is 14% less massive than the receding side and its rotation curve does 
not reach terminal velocity. The rotation curve of the receding (southeast) side however 
resembles that of a normal galaxy. 

Twelve high-latitude features were catalogued, of which four exhibit the signature of 
an expanding shell. There is some correlation of these features with features observed 
in the radio continuum from independent data. The most massive shell (Feature 10) 
extends out to 6h~ l kpc above the galactic disk. The radii and masses of these shells 
range from 1.2 h^ 1 — 2.0 hr 1 kpc and 1.6 x 10 7 — 5.5 x 10 7 h~ 2 Mq, respectively. We have 
investigated the possibility that the supershells could have been produced by external 
impacting clouds, but conclude that this scenario is unattractive, given the age of the 
shells, the isolation of the galaxy, and the lack of any observed features sufficiently 
massive to form the shells in the vicinity of the galaxy. Therefore, an internal origin is 
suggested. Since the implied input energies from supernovae are extremely high (e.g. 
from 1.4 x 10 53 /i~ 2 — 7.4 x 10 55 /i~ 2 ergs, corresponding to 400 — 74,000 supernovae), 
we suggest that some additional energy (e.g. from magnetic fields) may be needed to 
produce the observed supershells. 

Subject headings: 
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1. Introduction 



More than a decade ago, Heiles found HI shells, 
supershells and "worms" in the Milky Way Galaxy 
(Heiles 1979, 1984). These structures are seen pro- 
truding from the plane of the Galaxy in the z-direction. 
Diameters of the shells and supershells range from a 
few tens of pc to a few kpc. In some cases, the shells' 
diameters change with velocity, indicative of expan- 
sion. Similar features are also found in some exter- 
nal galaxies (e.g., NGC 5775, Irwin 1994; NGC 4631, 
Rand and van der Hulst, 1993) either as shells, par- 
tial shells or extensions. M31 is also observed to have 
hundreds of HI "holes" , regions devoid of neutral hy- 
drogen, which are possibly cavities within shells or 
supershells (Brinks and Bajaja 1986). Several other 
example also exist, but the list is still short. 

Two possible explanations have been proposed: 1) 
the effect of stellar winds and supernova explosions 
and 2) collisions of clouds with a galactic disk (see 
Tcnorio-Tagle and Bodcnheimcr 1988 and references 
therein). Considerable effort has been expended in 
the first category beginning with the "galactic foun- 
tain" model of Shapiro and Field (1976). Norman 
and Ikeuchi (1989) then extended this to develop the 
"Chimney Model" in which they propose that high- 
latitude features are formed due to clustered super- 
nova explosions. The hot gas in the disk flows to 
the halo via "chimneys" with walls of HI gas, cools 
at high-latitude and subsequently falls back onto the 
disk. This way, a circulation of mass, energy and 
momentum is set up between the disk and the halo. 
It has been argued ( Rand and van der Hulst 1993 
and others) that kpc scale features observed in our 
own and other galaxies are too large and require too 
many input supernovae (e.g., >10,000 in some cases) 
to be explained in this way. Thus there exists an 
energy problem for large HI supershells if they are 
created by supernovae and stellar winds. The cloud 
collision model, in part a response to this problem, 
on the other hand, requires the presence of a com- 
panion. Thus the discovery of HI shells in isolated 
systems would strongly favour an internal origin. We 
have therefore targeted one such galaxy, NGC 3044, 
which we felt to be a good candidate for HI extensions 
(see below) and has no nearby companions. 

NGC 3044 (Fig. 0, from the ST Scl Digital Sky 
Survey Q) is an edge-on, SBc galaxy in the Leo Cloud 
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Fig. 1. — Optical image of NGC 3044 obtained from 
the ST Scl Digital Sky Survey. 



( [TuUy 1988j ). Solomon and Sage (1988) classify it as a 
Type I galaxy (i.e, second lowest level of interaction), 
which has companions within IOD25 (0.23/i _1 Mpc 
in projection, where h is given by i/o/lOO), on the 
Palomar Sky Survey (PSS) print, with a velocity dif- 
ference of less than 1000 km s _1 and showing no mor- 
phological disturbance. A dust lane (a = 9 ft 51 m 2!0, 
S = 1°49'23'.'5) appears prominently on the north- 
western (NW) side of the galaxy, obscuring part of the 
disk. Above the dust lane, a faint feature protrudes 
at an angle away from the major axis appears to be 
the source of the classification by Sanchez-Saavedra 
et al. (1990) of the galaxy as having a "barely per- 
ceptible" optical warp pointing in a counter-clockwise 
direction (inverse integral sign) . On the southeastern 
(SE) side, the galaxy's optical disk fades out more 
gradually than the NW end, suggesting a possible 
asymmetric distribution of matter. 

NGC 3044 was originally observed as part of a sur- 
vey to search for galaxies with extended radio con- 
tinuum halos. The results (Sorathia 1994) show the 
galaxy to have extended radio continuum emission 
to distances as far as 8 kpc from the plane. Since 



this galaxy is also infrared bright (Soifcr et al. 1987), 
it therefore appeared to be a good candidate for 
finding high-latitude neutral hydrogen arcs and fil- 
aments. Throughout this paper, we use a distance 



lr The Digitized Sky Surveys were produced at the Space Tele- 



scope Science Institute under U.S. Government grant NAG W- 
2166. The images of these surveys are based on photographic 
data obtained using the Oschin Schmidt Telescope on Palomar 
Mountain and the UK Schmidt Telescope. The plates were 
processed into the present compressed digital form with the 
permission of these institutions. 
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to NGC 3044 of 16. 14/i -1 Mpc (see §fj). At this 
distance, 1" = 0.078/i -1 kpc. 

In the following, we present the first high resolu- 
tion HI observations of NGC 3044. Thus, aside from 
searching for high latitude HI, we also model the HI 
distribution in the disk, deriving global parameters. 
Observations and data reduction are given in the 
HI distribution, velocity field and the global HI pro- 
file of the galaxy in §|3|. §|] provides the details of 
modelling the galaxy as well as a discussion of the 
modelling results. The observed asymmetry and the 
energetics of the HI supershells are given in §||. Fi- 
nally, a summary can be found in §pl 

2. Observations and Data Reduction 

NGC 3044 was observed with the Very Large Array 
(VLA)0 on June 19-20, 1993 in the C configuration. 
The primary flux calibrators were 3C48 and 3C286, 
having flux densities of 16.09 Jy and 14.88 Jy respec- 
tively. The secondary calibrator, 0922+005 (flux den- 
sity 0.75 Jy), was observed about every half hour. The 
data were Hanning-smoothed on-line. Bandpass cali- 
bration was done using the primary flux calibrators. 
33 line-free channels were averaged together and sub- 
tracted from all 63 channels to obtain spectral line 
data alone. The NRAO AIPS (Astronomical Image 
Processing Software) routine "MX" was used for map- 
ping and deconvolution of the data cube (using the 
"Cleaning" algorithm). Two cubes were obtained, 
one using natural weighting and the other using uni- 
form weighting. Correction for primary beam attenu- 
ation was applied to both cubes. Parameters pertain- 
ing to the observations and mapping of NGC 3044 are 
listed in Table @. 

3. Results 

3.1. The HI Distribution 

3.1.1. The HI Channel Maps 

The distribution of neutral Hydrogen in veloc- 
ity space can most easily be discerned from channel 
maps. The naturally weighted channel maps in Fig- 
ure ||(a) show clearly that the SE side of the galaxy is 
receding and the NW side is approaching. Since the 
dust lane occurs most prominently along the south 
edge of the galaxy, this must be the closest edge to 



us. Therefore, any trailing spiral structure observed 
should have an inverted "s" shape. 

There is distinct evidence of HI arcs and exten- 
sions and high latitude features away from the plane 
of the galaxy. For example, on the SE side, a dis- 
connected feature at a = 9 h 51 m 7r7, 8 = l°47'21'/0 
reaches a height of 8.4/i -1 kpc in the 1422 km s _1 
channel and an impressive high-latitude extension is 
found on the NW side at a = 9 fc 51 m 3!3, 5 = l^O'O'.'O 
from 1110 km s -1 to 1235 km s _1 . It reaches a pro- 
jected height of 7.1ft. -1 kpc from the midplane. These 
features and others will be discussed in more detail in 



15.2. 



The channel maps also reveal the asymmetric dis- 
tribution of gas in the galaxy, which we highlight in 
Figure |^(b) , by superimposing selected channel maps 
of equally red- and blue-shifted velocities with respect 
to systemic. The HI distribution of the receding side 
(dotted contours) is much more elongated than that 
of the approaching side (solid contours) and the ap- 
proaching side is more "active" displaying the largest 
high-latitude extensions above the HI disk. The disk 
appears to be fairly straight at high velocities, al- 
though the outer contours (1.5 - 3<r) on the NW side 
dip below the major axis, probably due to an exten- 



sion. The asymmetry will be discussed in £5.1 



2 Operated by Associated Universities, Inc. under contract with 
the National Science Foundation. 



3.1.2. The Column Density Maps 

In Figure |^(a), the velocity-integrated column den- 
sity (moment zero) map of NGC 3044 for the uni- 
formly weighted data is superimposed on the optical 
grey-scale image from the Digitized Sky Survey. The 
naturally weighted column density map (contours and 
grey-scale) are shown in Figure ||(b) . In Figure ||(a) , 
one extension is p art icularly obvious - labelled F10 
(feature 10, see §|5.2j). The disk's column density 
peaks at a = 9 h 5l m T.7, S = 1°48'48'.T. Besides this 
peak, there are three other peaks along the major 
axis. The two outermost ones are roughly symmetri- 
cally located with respect to the central peak, both at 
a distance of about 10/i -1 kpc (126") on cither side. 
The NW peak is 2.5 times stronger in intensity than 
its counterpart on the SE side, being 8.5 x 10 21 cm~ 2 
and 3.3 x 10 21 cm -2 , respectively. The peak of the HI 
distribution does not coincide with the peak of the op- 
tical image but is offset by about 43" along the major 
axis to the SE. The optical distribution extends to the 
outer SE HI peak, but with very low intensity. Inter- 
estingly, there is another HI column density enhance- 
ment (Njji — 7.6 x 10 21 cm -2 ) between the central 
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Table 1 

HI Observation and Map Parameters 



Parameter Value 



VLA configuration C 

Observing date 1993 June 

On-source observing time 8.75 hr 

Approximate largest scale visible a 420" 

Primary beam FWHM 31.'5 

Band center b 1318 km s _1 

Total bandwidth 1331 km s _1 

Channel width (=resolution) 20.8 km s _1 

Synthesized beam parameters: 

Uniform weighting 13'.'6x 13'.' 5 @ PA=-78?19 

Natural weighting 20'.'7x20'.'2 @ PA=-29?54 

Root mean square map noise: 

Uniform weighting 0.61 mJy beam^ 1 

Natural weighting 0.45 mJy beam^ 1 

Rayleigh- Jean conversion factor (Tf, / S) : 

Uniform weighting 3.30 K(mJy beam -1 ) -1 

Natural weighting 1.45 K(mJy beam -1 ) -1 



a From Perley, 1994 
b Heliocentric, optical definition. 
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Fig. 2. — (a) Naturally weighted velocity channel maps of NGC 3044. Contour levels are at 1.0(1. 5c), 1.3, 1.9, 3.2, 
6.4, 12.8, 19.2, 25.6, 32.0, 38.4, 44.8 mJy/beam. The channel width is 20 km s" 1 centred at the velocity which 
appears at the upper left or right corner of each frame. The synthesized beam is shown at the lower left corner of 
the first frame. 



peak and the outer NW peak, directly below feature 
10. This smaller peak seems to coincide with a region 
of low surface brightness on the optical image, which 
could be due to absorption by a dust lane. 

Measured from the outermost contour of the nat- 
urally weighted map [Fig. [§(b)], the HI disk spans a 
length of 5.'7 or 2.3xR 2 5 along the major axis. Fea- 
ture 10 is clearly visible and more extended than in 
Figure ||(a) along with various other extensions. F10 
extends out to about 6ft." 1 kpc above the plane of the 
galaxy in Figure ||(b). In fact, Figure ||(b) shows an 
extensive disturbance in the region near feature 10. 
We show below (see § |5.2| ) that feature 10 is indeed 
an expanding feature. Other protruding features are 
also visible in Figure ||(b). For example, the hole at 
a = 9 fe 51 m 4!0, 6 = 1°48T0" and numerous protru- 
sions along the northern edge of the galaxy can be 
seen. 

3.2. The Velocity Field 

The naturally weighted velocity field of NGC 3044 
is shown in Figure || The kinematic major axis is 
parallel to the optical major axis to within one de- 



gree out to a radius of about 50" (4/i _1 kpc). The 
kinematic minor axis is not perpendicular to the in- 
ner kinematic major axis, this is typically an indica- 
tion of a bar structure (Bosma, 1981), which is con- 
sistent with the optical classification of this galaxy. 
Beyond about 1.'7 (8/i _1 kpc) on the NW side, the 
velocity field becomes disturbed. For example, there 
is a closed 1160 km s _1 contour at the position of 
F10. On both sides, beyond Ah" 1 kpc, the kinematic 
major axis bends slightly towards the south. There- 
fore the major axis does not resemble the more com- 
monly seen kinematic warp in other edge-on galaxies 
and we interpret the disturbance near feature 10 as 
high-latitude features instead of a simple warp. 

In Figure 4, we plot the systemic velocity of the 
galaxy, obtained from the global profile, on the ma- 
jor axis of the galaxy (represented by Q at a = 
9^51 m 5!9, 6 = 1°48'58'.'7). This position is signifi- 
cantly different from the position of both the column 
density peak (labelled o in Fig. 4) as well as the opti- 
cal centre (labelled +), so that the kinematic centre 
is between the column density peak and optical cen- 
tre. However, the kinematic centre does coincide with 
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Fig. 2. — (b) Selected velocity channels symmetrical with respect to the systemic velocity are superimposed. The 
rotation curve peaks at 150 km s _1 (Table [|). Dotted contours are positive in velocity and solid contours are 
negative. The cross represents the kinematic centre and the major axis is shown as the horizontal line across. 
Contours are as in Fig. 2(a). Note that the images have been rotated so that the major axis appears horizontal. 



the peak of the radio continuum map (represented by 
A). The kinematic centre given by our model (la- 
belled □ ) is also shown. This will be discussed more 
fully in §4.1. Here, we identify the radio continuum 
ccni trc aa the 'true' centre of the 



galaxy. 



] Notice all velocity contours have bends and kinks 
in them, especially near the outer edges of the con- 
tours. Such features are usually explained by shock 
fronts occurring at the inside edge of spiral arms 
as a density-wave perturbs the local velocity field 
(Visscr 1980). However, due to the high inclination of 



NGC 3044, this is difficult to verify from the column 
density map. At the location of feature 10, the veloc- 
ity field is inconsistent with normal galactic rotation; 
instead it is more blueshifted. If feature 10 is actu- 
ally located at a larger distance (so that its distance 
from the nucleus in Figure ^ is just a foreshortened 
distance in projection), then we would expect its ra- 
dial velocity to be redshifted. It is clear that feature 
10 does not follow the general flow of material in the 
disk below it. 



3.3. The HI Global Profile 

Figure | shows the global profile of NGC 3044 
from the naturally weighted cube. The profile agrees 
within error with previously published data (Staveley 



Smith and Davies 1988, 12' beam) 



The global properties of NGC 3044 are tabulated 
in Table ||. Row 1 gives the systemic velocity from the 
midpoint of the 20% intensity level. Row 2 gives dis- 
tance deduced from the systemic velocity corrected 
to the reference frame of the 3K background radi- 
ation (se e the Third Reference Cat alogue of Bright 
Galaxies [ de Vaucouleurs et al. 1991 1, hereafter RC3). 
Row 3 gives the velocity width at the 20% level un- 
corrected for inclination due to the high inclination 
for this galaxy (i = 84° ± 2°, |Bottinclli et al. 1984[ ). 

The integrated flux density of NGC 3044 is given 
in row 4. Our value agrees with both Staveley-Smith 
and Davies (1988) and Krumm and Salpeter (1980). 
Row 5 shows the estimated HI mass of the galaxy 
using the equation M m = 2.35 x 10 5 D 2 f S ■ dV 
where Mhi is in Mq , D is in Mpc and J S ■ dV is 
in Jy-km s . This equation is valid under the as- 
sumption of optical thinness. Haynes and Giovanclli 
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Table 2 

Global HI Properties of NGC 3044 



Row Parameter Value 



(1) 


V sys (km s _1 ) 


1287 ± 10 


(2) 


D H ubbie {hT 1 Mpc) 
AV 20% (km s" 1 ) 


16.14 


(3) 


351 ± 10 


(4) 


JS-dV (Jy-kms) 
M H i(h~ 2 10 9 Mq) 


48.2 ± 3.6 


(5) 


3.0 ± 0.2 


(G) 


M hi /L b {M q /LqY 
Mrih^lO 11 Mq) b 


0.22 ± 0.02 


(7) 


1.1 ± 0.1 


(8) 


M T /L B (h Mq/L q ) 
Mm/Mrih- 1 ) 


8.1 ± 0.4 


(9) 


0.027 ± 0.002 



a L_B is calculated using the total "face-on" apparent blue mag- 
nitude (Bj) given in the RC3 and using a value of +5.48 for the 
absolute blue magnitude of the Sun. 

b Radius is taken to be the maximum extent of the third con- 
tour in the naturally weighted position-velocity diagram (Fig. ^) 
and V ro t is the average of the maximum blue- and red-shifted ve- 
locities with respect to the systemic velocity measured using the 
same contour. The third (4. 5c) contour is used so as to avoid the 
protruding feature on the approaching side of the p-v diagram (see 
Fig. which probably does not reflect the maximum blueshifted 
rotational velocity. 
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Fig. 3. — (a) Uniformly weighted column density map superimposed on the optical Digital Sky-Survey image 
(Fig. |l|). The map is Hanning smoothed in velocity space using 3 channels and spatially using a Gaussian of 
FWHM=20". A cutoff at the la level is applied. Contour levels are at 1.8, 3.0, 4.8, 6.0, 12.0, 30.1, 48.1, 60.0, 72.2, 
90.2, 108.3, 120.3 xlO 20 cm~ 2 . Feature 10 is labelled as F10. 



(1984) investigated the effect of HI self- absorption on 
HI integrated flux as a function of morphological type 
and galaxy inclination. They concluded that for Sc- 
type galaxies with inclinations like NGC 3044, the 
correction factor is about 1.3. Thus Mhi could be up 
to 30% higher than the value in Table ||. The peak 
brightness temperature for NGC 3044 is 83.1 K in the 
13" beam. Thus optical depth effects should not be 
strong. The HI mass of the receding side is 14% larger 
than the approaching side while the uncertainty due 
to the error bars in the global profile is only 8%. This 
asymmetry is obvious in Fig. ^| and will be discussed 



5.1 



Row 6 in Table | gives the HI mass to blue lumi- 
nosity ratio. Row 7 gives the estimated total mass of 
the galaxy within the outermost detectable HI radius. 
The usual equation assuming a spherical geometry is 
used [M T = 2.33 x 10 5 RV? ot (M Q )}. The quoted un- 
certainty reflects only the asymmetry of the receding 
and the approaching sides of the galaxy, which domi- 
nates the errors. Using a spherical geometry may re- 
sult in an overestimate of the true total mass within 



R if the total mass is dominated by the disk com- 
ponent. The overestimate will be no more than 40% 
(Lequeux 1983). Row 8 gives the ratio of total mass 
to blue luminosity and row 9 gives the fractional HI 
mass. 

The values of Table ^ were compared with those of 
Roberts and Haynes (1994) for the Local Supercluster 
(galaxies with v < 3000 km s -1 ), correcting to their 
value of H . With an additional slight correction for 
the method of calculating Mt, we find that all values 
for NGC 3044 are typical for its Hubble type. 

4. Data Cube Modelling and Results 

For edge-on galaxies, the usual first and second 



moment analysis (van Gorkom and Ekers 1994) can- 
not be used successfully to obtain density and velocity 
distributions because the flux density at each pixel re- 
sults from an integration along the line of sight which 
spans many galactocentric radii. 

We instead use the technique of Irwin and Seaquist 
(1991) and Irwin (1994) which models the intensity at 
every pixel given the radial and perpendicular den- 
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Fig. 3. — (b) As in (a) but using naturally weighted cube. Column density contours are superimposed on grey-scale 
map. Contour levels are at 1.6, 2.6, 3.2, 4.0, 5.3, 13.2, 26.4, 39.7, 52.9, 66.1, 79.3, 92.5, 97.8xl0 20 cm- 2 . 



sity distributions and the form of the rotation curve. 
These curves are parameterized, and the parameters 
are fit via a non-linear least squares algorithm. The 
parameters are a subset of RA and DEC of the nu- 
cleus, the position angle (PA) and inclination (j) of 
the galaxy, the systemic velocity (V sys ), the galacto- 
centric radius (R max ) at which maximum rotational 
velocity (V max ) occurs, a rotation curve shape in- 
dex (m), the peak volume density of the distribution 
{n m ax), the density scale length for the radial dis- 
tribution (r Q ) and the scale height (z Q ) of the ver- 
tical distribution. The model intensities are then 
smoothed to the spatial resolution of the data and 
the residuals (data minus model) are found. The best 
result is considered to be the model producing the 
lowest residuals. 

The trial density distributions are exponential or 
Gaussian for both the radial and z-axis distributions 
centered at the galaxy centre and midplane, respec- 
tively. In addition, a Gaussian or exponential ring 
distribution in the plane of the galaxy centered at 
a radius of R Q was also tried. The choices for the 
form of the velocity curve are the Brandt curve and a 
user specified numerical rotation curve which we take 
from the velocity-position diagram along the major 
axis (Fig 6) (allowing the amplitude to vary). The 
model also allows for a velocity dispersion. Thus 



the geometrical parameters (particularly the inclina- 
tion which cannot normally be found without an as- 
sumption of the scale height), the velocity distribution 
and the density distribution (radial and vertical) can 
be found. For NGC 3044, we model the uniformly 
weighted data cube as it has a higher spatial resolu- 
tion than the naturally weighted cube, yet the inte- 
grated flux density (42±6 Jy-km/s) still agrees with 
that of the naturally weighted cube (see Table ||). 
There are 2248 independent data points of real pos- 
itive emission, hence the free parameters should be 
well constrained. 

The best fit results are shown in Table |^ for the 
entire galaxy (column 2) and the receding (column 3) 
and advancing sides (column 4). Column 1 is a list 
of free parameters. Brandt rotation curves are used 
for all the models in this table, since they resulted 
in better fits in all trials than a user-specified numer- 
ical curve. We shall hereafter use notations such as 
'RG(GS)-EX' to represent a model having a Gaussian 
ring distribution in the plane and an exponential dis- 
tribution in the z-direction. The uncertainties quoted 
are either standard deviations of the results of all trial 
models (including those not shown in Table |^) or for- 
mal lcr errors from the fit of the best model, whichever 
is larger. The following subsections discuss the model 
results. 
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Table 3 
Model Parameters 



Parameter 


Whole Galaxy 


Receding Side 


Approaching Side 




RG(GS)-GS 


GS-GS 


RG(GS)-EX 




ring @ 30" 




ring @ 104" 


(1) 


(2) 


(3) 


(4) 


a o (1950) 


9 fe 51 m 6?15 ± 0!07 






(5 o (1950) 


l°48'51'/2 ± Of! 4 






PA(°) 


112.95±0.10 


113.45±0.04 


111.4±0.3 




84.9±0.4 


85.3±0.3 


86.1±0.2 


V sys (km s 1 ) 


1298±4 






V max (km s _1 ) 


150±3 


149±2 


151±3 


Rmaxih^ 1 kpc) 


6.4±0.4 


4.9±0.2 


9±4 


Roih' 1 kpc) 


2.3 




8.2 


m 


0.99±0.18 


0.64±0.15 


0.48±0.26 


n m ax(h cm- 3 ) 


0.38±0.15 


0.40±0.18 


0.33±0.16 


r (/i _1 fcpc) 
o (/i _1 kpc) 


5.3±0.7 


6.4±0.8 


1.9±0.2 


0.42±0.04 


0.49±0.04 


0.45±0.03 


r 0i i(/i _1 kpcf 
u v (km s _1 ) 


1.20±0.02 




2700±29 


24.3±2.2 


19.9±2.2 


17.7±2.2 



a Inner Scale Length for r < R a . 
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Fig. 4. — The naturally weighted velocity field of 
NGC 3044. Contours are 20 km s _1 apart starting 
from 1130 km s _1 to 1410 km s" 1 . Also included are 
contours of 1160 and 1420 km s _1 . The kinematic 
major axis is indicated by the straight line across. 
Symbols represent the different peaks and centres de- 
scribed in the text. + is the optical centre, A is the 
radio continuum peak, o is the column density peak, 
O is the centre based on the global HI profile and □ 
is the model centre. 

4.1. Galaxy's Geometry 

The position of the nucleus of the galaxy is well 
constrained amongst the different models giving a 
standard deviation of 0!07 in RA and 0'.'4 in DEC. 
This is not surprising since the kinematic centre of 
a galaxy does not depend on the form of the density 
distributions either parallel or perpendicular to the 
plane. As mentioned before, the column density peak 
is shifted towards the SE (receding side) (§f^2j). As 
part of the modelling process, we forced the model's 
nucleus to coincide with the column density peak and 
the radio continuum peak separately. These resulted 
in worse models for all combinations of density dis- 
tributions implying real offsets of these peaks from 
the model's centre. The kinematic centre as given by 
the best model is 2'.'3 east and 8'.'8 south of the ra- 
dio continuum peak (^Sl^O, l°49'0'/0, see Fig. |). 
As the receding side is more massive, it is natural 
that the model's centre is being weighted towards this 
side, hence the offset of the model's centre from the 
radio continuum centre. The peak of the column den- 
sity map is 23'.'3 east and 3'.'1 south of the modelled 
galaxy centre. These differences in the peak positions 
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Fig. 5. — The naturally weighted global profile for 
NGC 3044. Error bars are at the la level. 

are significant. The P-V diagram of the major axis of 
the galaxy (Fig. ||) also shows an asymmetry in the 
sense that the receding side of the galaxy extends fur- 
ther. In addition, the flat part of the rotation curve 
shows up more prominently on this side. It there- 
fore appears that either the galaxy's HI emission has 
been 'truncated' on the approaching side or the re- 
ceding side emission has been 'stretched out'. We feel 
that the former is probably true as the receding side 
of the P-V diagram more closely resembles the P-V 
diagrams of other giant spirals. 

The position angle and inclination of the galaxy are 
also well determined. Various models resulted in stan- 
dard deviations of 0?10 and 0?4 for these two parame- 
ters, respectively. The position angle agrees with that 
given in the RC3 (note that there is a mistake in the 
RC3, the position angle listed should be 113° instead 
of 13°). The inclination from the model (84?9) agrees 
with Bottinelli et al. (1984) and Staveley-Smith and 
Davies (1988) who obtain i=84±2° and 90±13° re- 
spectively. The former assumes an intrinsic oblate- 
ness of 0.15 while the latter uses 0.20. Our value is 
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Fig. 6. — Position- velocity diagram along the kine- 
matic major axis of NGC 3044. The data cube used 
is naturally weighted. Contour levels are 0.8(1.8cr), 
1, 2, 5, 10, 15, 20, 25, 30, 35, 40 mJy/beam. The 
cross indicates the radio continuum centre placed at 
V sys based on the global HI profile. Note that at 
the galaxy's declination and P. A., I s corresponds to 
16'.'3 along the major-axis. Same for all similar figures 
hereafter. 



independent of the intrinsic axial ratio, and should be 
an improvement over the previous estimates. 

4.2. The Rotation Curve 

The systemic velocity given by the best model 
agrees with that measured from the global profile of 
the data (Table ||) to within errors. Fixing the sys- 
temic velocity at a lower value resulted in worse mod- 
els for all trials. The rotation curve parameters, V max 
and Rmax, are also well-defined, with the spread in 
Vmax within 10% of the velocity resolution while the 
spread in R max is determined to 40% of the HPBW. 
The shape of the rotation curve is determined by the 
Brandt index, m, and the spread in this parameter is 
only 18% amongst all trials. When a velocity disper- 
sion is included in the models, results for every model 
improved significantly. For the best model, the opti- 
mum FWHM of the Gaussian smoothing function is 



found to be 57.2 km s" 1 (2.75 channels) corresponding 
to a velocity dispersion of 24.3 km s -1 . This value re- 
flects all non-circular motions along each line of sight, 
including any velocity variations due to spiral arms 
and is a global average. This velocity dispersion is 
comparable to that found in NGC 5775 (Irwin, 1994) 
and is about 3 times the velocity dispersion found for 
the Milky Way galaxy, which only measures the rms 
cloud velocity (Spitzer, 1978). 

4.3. The Density Distributions 

The model parameters related to the density dis- 
tributions vary significantly amongst different mod- 
els. This is expected as the parameters are de- 
fined differently between some of the models. While 
there is no obvious central hole in the column den- 
sity map (Figs. 3), the best results of the modelling 
are found using a ring distribution rather than a dis- 
tribution which peaks at the centre. From our best 
model (Gaussian ring in-plane and Gaussian in z- 
direction), the volume density at the radius of the 
ring (2.3/j _1 kpc) is found to be 0.38/i cm -3 . This 
is very similar to the average midplane density of 
our Galaxy between 4 to 8 kpc which is w0.35 cto~ 3 
( Burton and Gordon 197"s| ). The outer radial scale 
length is 5.3/i _1 kpc, or 20% of the length of the HI 
disk as measured by the outermost contour of the col- 
umn density map [Fig. ||(b)]. As a comparison, the 
Galactic HI density stays roughly flat from about 4 
to 10 kpc and falls off beyond that. The inner scale 
length is the least constrained, its value depends sen- 
sitively on the position of the ring. The uncertainty 
associated with this parameter in Table || is there- 
fore taken to be the 1-a uncertainty of the fit instead 
of the standard deviation amongst the different mod- 
els. Our best model gives the inner scale length of 
1.20ft." 1 kpc, about half the radius of the ring. 

The vertical scale height given by the best model 
is 0.42/i -1 kpc. This value is obtained as a global pa- 
rameter, an average over the entire disk of the galaxy. 
Note that in specific places, e.g., feature 10, the HI 
extends to 5.7h~ 1 kpc in projection. In the Galaxy, 
it is known that the HI distribution between 4 to 
8 kpc consists of two components, a central layer with 
a FWHM of 0.1 kpc, and a low- intensity, high- 
temperature component with a FWHM of «0.5 kpc 
(Lockman, 1984). It is perhaps more helpful to com- 
pare the vertical scale height of NGC 3044 with an ex- 
ternal galaxy. NGC 891 is also an edge-on, IR-bright 
spiral which exhibits high-latitude HI features, van 
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Table 4 

High-Latitude Arcs and Extensions 



Feature Right Ascension Declination 


Velocities 


Number of 


Height from midplane 


h m s o 1 II 


km s _1 


Channels 


h^ 1 kpc 



1 


9 51 9.8 


1 49 5.0 


1401.2 to 1442.8 


3 


2.9 


2 


9 51 6.9 


1 49 21.0 


1359.6 to 1422.0 


4 


3.5 


3 


9 51 5.0 


1 48 17.0 


1234.8 to 1359.6 


7 


6.0 


4 


9 51 10.6 


1 48 13.0 


1338.8 to 1297.2 


3 




5 


9 51 7.7 


1 47 21.0 


1380.4 to 1463.6 


5 


8.4 


6 


9 51 9.5 


1 49 5.0 


1255.6 to 1318.0 


4 


5.0 


7 


9 51 4.7 


1 48 17.0 


1276.4 to 1318.0 


3 




8 


9 51 6.9 


1 49 29.0 


1193.3 to 1276.4 


5 


3.7 


9 


9 51 5.3 


1 50 1.0 


1193.3 to 1214.1 


2 


6.6 


10 


9 51 2.6 


1 50 0.0 


1110.2 to 1234.8 


7 


7.1 


11 


9 50 58.9 


1 50 2.0 


1151.7 to 1214.1 


4 


4.9 


12 


9 50 57.8 


1 49 0.0 


1131.0 to 1193.3 


4 


5.2 



der Kruit (1981) modelled the thickness of its HI layer 
and found that the FWHM of the z-distribution in- 
creases with galactocentric radii. The FWHM varies 
from 0.32 kpc at a radius of 4.2 kpc to 1.89 kpc at a 
radius of 20.8 kpc with an average value of 0.9 kpc. 
Therefore, NGC 3044 appears to have a moderately 
thick HI disk compare to the Galaxy yet not as ex- 
tensive as that in NGC 891. To check the validity 
of a thick disk, we modelled the galaxy again by fix- 
ing the vertical scale height to a low value (0.1 kpc). 
This resulted in a significantly worse model fit to the 
data. In addition, Sorathia (1994) finds a large radio 
continuum scale height of 1.8±0.5 kpc for the 20 cm 
C-array data (see Fig. 9). Thus NGC 3044 has both 
a moderately thick HI disk as well as a thick radio 
continuum disk. 

4.4. Modelling of Receding and Approaching 
Halves of the Galaxy 

As we have noted before, the HI distribution of 
NGC 3044 is not symmetrical in the sense that the HI 
is more extended on the SE (receding side) and less 
so on the NW (approaching side). Modelling the re- 
ceding and approaching side of the galaxy separately 
enforces this picture. Table || lists the best fit param- 



eters for the two halves of the galaxy in comparison 
with the parameters given by the galaxy as a whole. 
Column 3 and 4 were obtained by fixing the nucleus 
position and the systemic velocity given by the best 
model for the whole galaxy (Table ||, column 2). 

The model fits the receding side of the galaxy bet- 
ter which shows that this side of the galaxy is better 
described by a smooth spatial and velocity distribu- 
tion of gas. This result simply reiterates what we 
see in the major axis rotation curve from the data 
(Fig. §). Although V max and R max of both sides 
agree within error, the uncertainty in R max is much 
larger for the approaching side. The Brandt curve in- 
dexes show the same trend. The difficulty in pinning 
down these parameters is due to the fact that on the 
approaching side, the rotation curve never reaches a 
terminal velocity. For the density distribution param- 
eters, the two halves of the galaxy differ significantly. 
The receding side is best fitted by a GS-GS distribu- 
tion while the approaching side prefers a RG(GS)-EX 
distribution. In the latter case, the location of the 
ring is at the outer density peak and the inner scale 
length is extremely high, hence the HI distribution 
inside this peak is essentially flat. A consequence of 
the large ring radius is the small outer scale length, 
as expected. These results for both sides, separately, 
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are consistent with the residual emission observed in 
Fig. [?] (see next section). Finally, the approaching 
side has a smaller velocity dispersion than the reced- 
ing side, which can also be seen from the narrower (in 
velocity) emission on the approaching side in Fig. 6. 

4.5. Comparing Model With Data 
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RIGHT ASCENSION (B1950) 

Fig. 7. — Residual column density map obtained by 
subtracting model from data. Solid contours repre- 
sent emission above a smooth distribution and dotted 
contours represent too much model emission. Con- 
tour levels are at 36.1, 30.1, 24.1, 18.1, 12.0, 6.0, -6.0, 
-12.0, -18.1 and -24.1 xlO 20 cm" 2 . Feature 10 is la- 
belled as F10 on the map. The symbols have the same 
meaning as in Figure || 

Figure shows the column density map of the 
residual cube for the global model result. Negative 
(dotted) contours represent regions where the amount 
of gas was over-estimated by the model and posi- 
tive (solid) contours represent regions of enhanced 
emission in the galaxy above a smooth distribution. 
The peak column density in the residual map is 30% 
of the peak column density in the data [Fig. ||(a)], 
this shows that the model is able to reproduce 70% 
of the underlying smooth distribution of the galaxy. 
There remain, however, pixel-to-pixel variations be- 
tween model and data due to the uncvcnncss of the 
distribution, as can be seen in Fig. 0. For example, 
feature 10 appears as positive contours, as do the two 
outer column density peaks [see also Fig. 0(a)]. The 
galaxy's nucleus derived from the best model does not 
coincide with the central peak of the zeroth-moment 
map. As a result, the model under-estimated the col- 
umn density at the data peak. 

In terms of velocity, the model reproduced the ve- 
locity field of the galaxy well for the most part. This 
is confirmed by subtracting the first moment map of 
the model from that of the data (Fig. ||). The rms of 
this residual velocity field is 16 km s , smaller than 



the velocity resolution of the data. There are regions 
in the galaxy where the model's velocity cannot re- 
produce the data, a good example is the region just 
under feature 10. In this region, the typical velocity 
excess is about -50 km s" 1 (i.e., model over-estimated 
the velocity) which represents the highest departure 
of the model from the data. This is consistent with 
the observation made in §3.2. which shows that the 



velocity field near feature 10 is more blueshifted for its 
location than expected from normal galactic rotation. 
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Fig. 8. — Residual velocity field (contours) obtained 
by subtracting the velocity field of the model from 
that of the data. Contours are at -55, -50, -45, -40, - 
35, -30, -25, -20, -10, -5, 5, 10, 15 and 20 km s _1 . It is 
superimposed on the grey scale image of the velocity 
field of the data for comparison. Grey scale intensities 
are shown as a band at the top. The cross marks the 
radio continuum centre of the galaxy and F10 points 
to the location of feature 10. 



5. The Remarkable Features of NGC 3044 

In the following, we discuss two of the most re- 
markable features of NGC 3044, namely, the HI asym- 
metric distribution and the high-latitude structures. 
These features are readily understood for an interact- 
ing galaxy but remain puzzling in a seemingly isolated 
galaxy like NGC 3044. 

5.1. Asymmetry 

The asymmetric appearance of the optical disk has 
been described in §|l| while the HI asymmetry is shown 
in the HI global profile (Figure 5) and the rotation 
curve (Figure |gJ) and was discussed in § |3.l| . In Fig- 
ure ^, we compare the HI distribution with the 20-cm 
radio continuum map of similar beam size obtained 
from the completely independent data of Sorathia 
(1994). There are remarkable correlations between 
the 2 maps, both show extensions of the intensity 
contours towards the NW as well as enhancements 
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in intensity at the most north-westernly peak of the 
HI column density map [a = 9 h 51 m 0!0, S = 1°49'35"8 
cf. Fig. ||(b)]. All the above points to a global asym- 
metry which manifested itself in various components 
of this galaxy. 

As indicated in § [l], there is no morphological evi- 
dence for a disturbance in this galaxy (Solomon and 
Sage 1988). What, then, could have caused the asym- 
metric distribution in NGC 3044? We note here that 
HI asymmetries in field spirals are quite common. For 
example, at least 50% of 1400 field galaxies were clas- 
sified as asymmetrical by Richter and Sancisi (1994). 
However, the origin of the asymmetry still eludes us. 
Zaritsky (1995) suggests that past mergers may be the 
cause of the optical asymmetries in a sample of about 
30 galaxies. Could the asymmetry in NGC 3044 be 
similarly due to a past merger, perhaps with a smaller 
companion? 

Mergers between gas-rich disks and satellite galax- 
ies (i.e., minor mergers) have been studied extensively 
via simulations (e.g., Quinn et al., 1993, Mihos and 
Hcrnquist, 1994 and Hcrnquist and Mihos, 1995). In 
general, the effect on the parent galaxy for assimilat- 
ing a satellite one tenth its own mass are: 1. massive 
gas in-flow to the nucleus causing a brief (~ 10 8 years) 
starburst phase; 2. heating of the stellar disk (i.e., in- 
crease in vertical scale height by a factor of a few); 3. 
flaring and warping of the stellar disk. In addition, de- 
pending on the satellite's initial density, its core may 
or may not survive the tidal stripping to arrive at the 
nucleus of the parent. Therefore, an obvious signa- 
ture that a merger has occurred is the existence of 
a double nucleus. In addition, Zaritsky (1995) finds 
a possible correlation between HI asymmetries and 
star formation rates using the set of Hi-asymmetric 
galaxies from Rix and Zaritsky (1995). The higher 
than normal star formation rates are consistent with 
the starburst phase expected after a minor merger. 

We do not yet have the data which could search 
for all of these effects. Star formation rates can 
be investigated, though. NGC 3044 is indeed clas- 
sified as infrared-bright, but this is based on the 
IRAS 60/zm flux density (Soifer et al. 1987). The in- 
frared luminosity calculated from the 60/im and the 
lOO^m flux densities given by Soifer et al. (1989) 
is 4.8xl0 9 /i -2 Lq and the massive star formation 
rate, following Condon (1992), is given by SFR(M > 
5M Q ) = 9.1 x 10- n -^a = O.Uh- 2 M Q /yr. These 

values show that NGC 3044 is only mildly starburst- 



ing at present (even if h = 0.75 is adopted, which 
gives SFR = 0.78 MQ/yr). That is, the massive star 
formation rate in NGC 3044 is comparable to that of 
the Milky Way (cf. 0.3 to 0.5 MQ/yr for a supernova 
rate of 1 every 5-80 years) and lower than that of 
M 82 (2.2 M Q /yr) by a factor of 5. Thus if a minor 
merger has occurred in NGC 3044, it likely occurred 
of order ~ 10 8 years ago or more. A study of the 
CO content is underway to investigate the star for- 
mation efficiency, determine whether the CO is also 
asymmetric and to search for evidence of an increase 
in molecular gas density towards the nucleus (Lee and 
Irwin, in preparation). Further observations (e.g. op- 
tical or IR searches for a double nucleus or a close 
examination of the optical disk for flaring) would also 
be useful. 

5.2. High-Latitude Arcs and Extensions 

5.2.1. Evidence of Expanding Shells 

The naturally weighted channel maps in Figure ^|(a) 
display numerous low intensity, high-latitude arcs and 
extensions away from the galactic disk. These arcs 
and extensions in NGC 3044 resemble the so-called 
"Heiles Shells" in our Galaxy. Table ^ lists the posi- 
tions, velocity ranges and the highest z-extent at the 
1.5 a level of the more prominent extensions/holes 
seen in Figure ||(a). Only features which show up 
in more than one consecutive channel are included. 
Some of these have very complex appearances and 
are probably the results of blended features along the 
same line-of-sight (e.g., feature 3). These features ap- 
pear to be distributed randomly along the disk of the 
galaxy. 

Features 4 and 7 are clearly "holes" and are most 
obvious at 1339 km s _1 and 1297 km s -1 , respectively 
[see Fig. ||(a)]. Feature 6 straddles either side of the 
systemic velocity (1287 km s _1 ), hence we may be 
seeing both the receding and the approaching caps 
of an expanding shell on the northeast side of the 
galaxy (see below). Feature 10 is the most massive 
(see Table |J) extension in the list and spans 7 velocity 
channels. At 1172 and 1152 km s _1 , feature 10 bends 
towards the east at high latitude so that the feature 
looks like an arc. Feature 5 reaches 8. Ah^ 1 kpc above 
the midplane (most obvious at 1422 km s _1 where it 
could be disconnected from the disk), making it the 
most extended feature of all. 

If any of these features are actually expanding 
shells or partial shells, they should appear as rings 
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Fig. 9. — 20-cm C-array radio continuum contour map is superimposed on the natural weighted HI column density 
map. The radio continuum data has a beam size of 18'.'87 x 15'.'86 and contours are at -0.15, 0.15 (lcr), 0.30, 0.45, 
0.75, 1.05, 1.50, 2.25, 3.00, 4.50, 6.00, 7.50, 9.00, 12.00, 15.00, 18.00, 21.00, and 22.50 mJy/beam. A cross marks 
the radio continuum centre of the galaxy. Note that the HI disk has been "burnt-out" in order to show the faint 
high-latitude features. 



or partial rings in the position-velocity (P-V) slices 
across the centres of the features. Figure 10(a) and 
(b) show the zeroth-moment maps of features 4, 7, 
10 and 12 and their associated P-V diagrams (nat- 
ural weighting data). Each slice is an average of 2 
pixels (8") parallel to the major axis. In each P-V 
slice shown in Figs. 10, we find expansion signatures 
at the position and velocity corresponding to these 
features. In all cases, the shells are not complete but 
are "open" at the sides toward lower density, both in 
P-V and in RA-DEC space. This is consistent with 
the Chimney models for a blow out case (see below). 
It could also indicate that the shell has more readily 
fragmented in low density regions. 

All the features in Table [| reach a projected height 
of at least 5 times the HI scale-height derived from 
our model (see Table |. It is clear that, whatever the 
mechanism which created them, it must be energetic 
enough that the features are able to break through 
the moderately thick HI disk to reach the halo (this 
is termed blowout in Hcilcs 1990[ ) . For a blowout case, 
mass and momentum are injected into the halo. The 
total mass of all features measured is 2 x 10 8 h~ 2 Mq 
(a lower limit, see next section) or 7% of the total HI 
mass of the galaxy. If we assume that 10% of this mass 



eventually reaches the halo, and the typical lifetime 
for these features is about 3 x 10 7 /i _1 years (see r s h 
in Table |^ and the next section), then a lower limit 
to the mass injection rate for cold neutral gas, M, 
is 0.7/i _1 Mq/vt. This is in rough agreement with 
numerical calculations utilizing the ISM parameters 
of the Galaxy flHeiles 1990j). 



5.2.2. Super shell Parameters 

Table | lists the masses (M s h, column 2), expansion 
velocities ( V s h> column 3), radii {R s h, column 4), pro- 
jected galactocentric locations (D s /i, column 5), am- 
bient densities in the mid-plane at D s h (rio, column 
6) , and kinetic energies (£* = \M s } l V^ h , column 7) 
of the expanding shells. The masses are obtained by 
summing the fluxes associated with each features in 
the velocity channels listed in Table [|. These masses 
are lower limits since part of the shells are most cer- 
tainly embedded within the HI disk and are not in- 
cluded. For holes 4 and 7, the masses are not listed as 
the thickness of the shell cannot be determined with 
confidence. V s h is taken to be half the total veloc- 
ity range which again may underestimate the actual 
range. Radii of the shells are measured from the 1.5<r 
contours of the P-V diagrams. D s h is an average of 
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Fig. 10. — (a) (i)Zeroth-moment of the velocity channels showing features 4 and 7. Smoothing is done as in Figs. [| 
(ii) and (iii) show position-velocity (P-V) slices parallel to the major axis across feature 4 and 7, respectively. Each 
slice is an average of 2 pixels (8"). Contours are at 0.64 (ler), 1, 1.3, 1.6, 1.9, 3.2, 5.1, 6.4, 7.7 mJy/beam. The 
height above and below the midplane is shown in each P-V panel. 
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Fig. 10. — (b) as in (a) but for feature 10 and 12. 
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the measurements made from the channel maps and 
the P-V diagrams. We use the modelled density pro- 
file in the plane of the galaxy (see column 2 of Table ||) 
to estimate the ambient density at the galactocentric 
distance of a given shell. 

The kinematical age of a shell, assuming constant 
expansion velocity, is r a h = R s h/V s h- For the four 
shells listed in Table ||, we expect V s h to be underes- 
timated (see above) and therefore T s h may be overes- 
timated. The resulting kinematical ages range from 
24 million years for shell 10, to 52 million years for 
shell 7. These ages are comparable to the expected 
lifetime of an OB association. 

For an expanding shell which is formed from a one- 
time energy injection, such as from supernovae and is 
now in the radiating phase of its evolution, numerical 
analysis by Chevalier (1974) shows that the energy 
injected is given by 

E E = 5.3xlO^R^V^ 

where the variables have the same meaning as before 
(R jit is in pc and Van in km s -1 ). Values of Ee for tho 



exp anding shells in NGC 3044 are listed in column 8 
of Table |. As Table | shows, since the energies are 
all upwards of 10 53 ergs, it is evident that all four 
shells are "supershells" , as defined by Heiles (1979, 
i.e. with energies > 3 x 10 52 ergs). 

Since the input energy requirements for supershells 
are significant, it is worth considering the errors on 
Ee- As indicated previously, the measured V s h is 
a lower limit since complete shells are not observed. 
Increasing V s h by a factor of 2, for example, would in- 
crease Ee by a factor of 2.6. We have also minimized 
R s h through our choice of Hubble constant (see de- 
pendence on h in Table [sj> . Thus a change to h = 0.75 
increases Ee by a factor of 1.8. If the features are ac- 
tually located at larger galactocentric distances than 
the projected radius, D s h, then the ambient densi- 
ties, no, would be lower than the listed values. In 
the extreme case of a shell actually occurring at the 
outer radius of the HI disk (R — 13.4/i _1 kpc), then 
no = 0.01ft- cm -3 , decreasing Ee by a factor of 30 
for feature 10, which has the highest uq in Table |. 
However, note that for Feature 10, = 3 x 10 54 
ergs. This value gives only the kinetic energy of ex- 
pansion of the feature and thus the input energy must 
be much higher (i.e. by at least a factor of 10 for an 
efficiency of 10%). Therefore, Ee for this feature is 
likely to be within a factor of 2 of the value listed 
in Table ||. This then implies that Feature 10 is at 



a position much closer to its projected galactocentric 
distance than near the edge of the disk. Given the ir- 
regular density distribution in the galaxy (§ [T^, Fig. 
7) , there are also possible variations in no for individ- 
ual shells. For example, Feature 7, with a projected 
center at the center of the galaxy, is the most extreme 
example since it could have an energy more than an 
order of magnitude larger than the table || value if no 
is determined by the fits to the two sides of the galaxy 
separately, rather than to the galaxy as a whole (Ta- 
ble ^). Overall, however, the irregularities in density 
tend to be within ~ 30% of the modelled density, in- 
troducing similar sized errors into Ee- Most of the 
assumptions we have made effectively minimize the 
values of Ee in Table ||. We conclude that the energy 
estimates, under the assumption of instantaneous in- 
put, are order of magnitude values. 

If supernova explosions are the source of the en- 
ergy, then between ss 400 to 70, 000 supernovae are re- 
quired to produce the shells in Table 0. Heiles (1979) 
found energies an order of magnitude lower for Galac- 
tic sh ells while shells in NGC 3079 (Irwin and Seaquist 



1990| ) and NGC 4631 ( |Rand and van der Hulst 1993[ ) 
require comparable energies 

(10^4 _ iq55 erg ) w h en 
scaled to h = 1. Note that NGC 3044 probably has 
many lower energy shells like those in the Galaxy and 
the fact that we detect only the supershells is merely 
a selection effect. Such a high energy requirement for 
the supershells is difficult to reconcile with input ener- 
gies typical of Galactic OB associations which contain 
only a few tens of stars of spectral types B0 and ear- 
lier. Ha luminosities of bright HII regions in some 
external late-type galaxies do suggest the existence 
of superclusters which contain thousands of super- 
novae (Heiles 1990) and therefore Ha observations of 
NGC 3044 may be helpful in this regard. 

There is certainly some dependence of the required 
input energy on adopted model. For example, if a 
slower, continuous energy injection is assumed over a 
typical cluster (« shell) age (several x 10 7 years), then 
we find energies which agree to within a factor of 4 
with the tabulated Ee for Features 4 and 10, but are 
at least an order of magnitude lower for the low den- 
sity Features 7 and 12 (cf. Vader and Chaboyer 1995). 
Also, since the supershells have achieved blowout, an 
understanding of these features would also benefit 
from numerical hydrodynamical modelling for such 
conditions (cf. Mac Low, McCray, & Norman 1989). 
Nevertheless, one cannot escape the fact that some of 
the observed high latitude features in NGC 3044 are 
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Table 5 

Energy Requirements of Expanding Supershells 



Feature 


M sh 


V sh 


Rsh 


D sh 


n Q 


E k 


Ee 




h~ 2 M Q 


km s _1 


hr 1 kpc 


h~ x kpc 


h cm~ 3 


h~ 2 ergs 


h~ 2 ergs 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



4 
7 
10 
12 



5.5xl0 7 
1.6xl0 7 



31.2 
31.2 
72.8 
41.6 



1.2 
1.7 
2.0 
1.5 



6.7 
0.0 
6.3 
10.1 



0.18 
0.01 
0.21 
0.04 



4.1xl0 54 
4.1xl0 53 
2.9xl0 54 7.4xl0 55 
2.7xl0 53 1.9xl0 54 



supershells which require large input energies to pro- 
duce. The kinetic energies of expansion, Ek, which 
are model-independent and essentially assume 100% 
input efficiency, are of order 10 53 - 10 54 ergs and these 
kinetic energies are underestimated for reasons out- 
lined above. For a more realistic situation in which 
the efficiency is 10% or less, the input energies must 
be at least an order of magnitude larger. 

5.2.3. Origin of the Supershells 

We have so far implicitly assumed that the shells 
are formed internally through the collective effects of 
clustered supernovae. Since the energy requirement 
are very high, however, the leading alternative sce- 
nario, i.e. that the supershells are produced from im- 
pacting external clouds, should also be considered. 

Impacting cloud models are an attractive way to 
explain large, energetic supershells because the re- 
sulting energies are a function of the infalling mass for 
which there are no hard limits in the case of galaxy- 
galaxy interactions. If we assume that the supershells 
in NGC 3044 are due to impacting clouds, however, 
we must consider where such clouds would originate. 
NGC 3044 has no nearby companion, nor is there 
strong evidence (apart from the asymmetry) for a pre- 
vious interaction. Moreover, a previous interaction, if 
it occurred, more likely occurred over timescales of 
order 10 s yr (a typical interaction timescale, see also 
§5.1). However, we find kinematical ages for the shells 
of a few x 10 7 yrs. Therefore, if the shells are pro- 
duced by infalling clouds, the situation is more likely 
one in which high velocity clouds (HVCs) exist around 



the galaxy and are continuously "raining down" . This 
also suggests that there should be evidence for such 
clouds around the galaxy now. 

One of the two supershell s found in NGC 4631 
(Rand and van der Hulst 1992) was modelled by Rand 
and Stone (1996) as a HVC impact structure using 
a 3-D hydrodynamical simulation. They found the 
most likely HVC that formed the supershell has a ra- 
dius of 500 pc and an HI mass of 1.2 x 10 7 Mq. Since 
NGC 3044 and NGC 4631 have remarkably similar 
kinematics, global HI distribution and shell param- 



eters (see Rand 1994), we can assume that infalling 
clouds of similar size and mass are required to form 
the supershells in NGC 3044. Converting these pa- 
rameters to a column density, we find (within one 
beam) a value of Njji = 5.2 x 10 20 cm -2 , which cor- 
responds to the fifth contour in Figure 3(b). There 
is clearly no evidence at the present time for such 
massive clouds in the vicinity of NGC 3044. Thus, 
the impacting cloud model may be reasonable for in- 
teracting galaxies like NGC 4631, but is much less 
attractive for an isolated system like NGC 3044. 

From Figure 9, there appears to be a correlation 
between high latitude radio continuum and high lat- 
itude HI features in the case of features 3, 6, and 
10. Indeed there appears to be an HI feature wher- 
ever a radio continuum feature exists. However, the 
converse is not true, as there are a number of HI fea- 
tures which do not (at the sensitivity limit of the ob- 
servations) have corresponding radio continuum fea- 
tures. The observations are not inconsistent, at least 
qualitatively, with the Chimney model ( Norman and| 
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Ikeuchi 1989). In this model, the HI features are walls 
of ambient material swept up by spatially correlated 
supernova explosions. Relativistic electrons then fun- 
nel through these "Chimneys" to reach the halo. It is 
difficult to tell, at this resolution and with the blended 
features, whether the radio continuum emission is in- 
terior to a specific chimney, or whether the emission 
is coincident with the walls of a chimney 

As the above analysis shows, the formation of su- 
pershells in NGC 3044 is far from clear. Both super- 
novae models and the cloud-galaxy collision model 
have difficulties. However, given the timescales for 
shell formation, the absence of companion galaxies 
near NGC 3044, and the lack of evidence for sur- 
rounding massive clouds, and the apparent correla- 
tion between high latitude HI and radio continuum 
features, we favour an internal origin for the super- 
shells. Several recent studies have shown, for exam- 
ple, that when magnetic fields are included (not just 
as a secondary effect) that large shells and blow-out 
can occur with more modest energy requirements (e.g. 
Kamayaet al. 1996; Frei et al. 1997). Therefore, more 
sophisticated models including effects like this may 
reconcile the energy requirements of supershells with 
fewer numbers of supernovae. 

6. Summary 

What has been presented in the preceding sections 
represents the first comprehensive HI study of the iso- 
lated, edge-on, infrared-bright galaxy NGC 3044. The 
IR brightness is partly due to proximity, since the 
massive star formation rate is a factor of 5 lower than 
that of M 82 (§5.1). 

The HI distribution in the galaxy was modelled us- 
ing all available data points, i.e. HI spectral lines were 
modelled over the entire galaxy in order to determine 
the best global density and velocity distributions in 
the galaxy as well as the best geometrical parame- 
ters. The resulting quantities are listed in Table ^. 
A Gaussian ring density distribution in the plane was 
found to be the best fit overall. An important result 
derived from the model is that the galaxy has a mod- 
erately thick HI disk with a Gaussian vertical scale 
height of 420/i~ 1 pc. This is consistent with previous 
observations which have also revealed a thick radio 
continuum disk. 

NGC 3044 shows a very asymmetric distribution 
optically, in HI, and in the radio continuum, with ex- 
cellent agreement between the radio continuum and 



HI column density distributions (see §5.1). In HI, the 
asymmetry is particularly apparent from the global 
profile, which shows a lopsided double-horned struc- 
ture. Otherwise, however, derived global parameters 
tell the story of a normal SBc galaxy. The column 
density map confirms the asymmetric distribution of 
the HI in the galaxy. This is clear from the 27" off- 
set of the HI peak to the east of the radio continuum 
peak. The major-axis rotation curve also shows the 
asymmetry in that the approaching side of the galaxy 
appears to be "truncated" before reaching terminal 
velocity while the receding side stays flat for about 
lhr x kpc after reaching terminal velocity. While the 
galaxy is asymmetric in HI, it also shows an extremely 
straight disk (i.e., no obvious disk warp). 

It is not clear what could have caused the asym- 
metry in this apparently isolated galaxy. It is likely, 
though, that the mechanism would have to act on the 
galaxy after its formation, since it is difficult to ex- 
plain how the galaxy would form asymmetrically in 
the first place. It is noted here that asymmetries in 
isolated galaxies are not uncommon, which indicates 
that a common mechanism such as past minor merg- 
ers may be responsible for it. If this is the case for 
NGC 3044, the merger likely took place over 10 s yrs 
ago. 

The velocity channel maps reveal a host of high- 
latitude HI extensions. We have catalogued a total of 
12 such features that appear in at least two consecu- 
tive velocity channels. These features are distributed 
uniformly across the disk of the galaxy and above and 
below the plane. We derive a lower limit for mass 
flow to the halo of M, = 0.7/i _1 M Q /yr. There is 
also some correlation between high-latitude HI fea- 
tures and high-latitude radio continuum features (§ 
5.2.3). 

Four of the high latitude features exhibit the sig- 
nature of an expanding shell and have radii of order a 
few kpc. Based on the assumption that the shells were 
formed by supernovae energy input, the required in- 
put energies are in the range 4 x 10 53 to 7 x 10 54 h~ 2 
ergs. These values vary somewhat depending on the 
adopted model, but it is clear, even from the kinetic 
energies of expansion of the individual supershells, 
that large numbers (e.g. up to tens of thousands) 
of correlated supernovae are required to form them. 
This means that NGC 3044 would have to harbour 
a number of such superclusters (such as the cluster 
R136 in the 30 Doradus region in the LMC), some- 
thing which is not supported by the current massive 
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star formation rate. 

The alternative model, that of impacting external 
clouds, is even less attractive as a formation mech- 
anism for the supcrshclls because of the absence of 
companions, the fact that the supershells are young 
(of order the age of a typical OB association) and 
because HI clouds which are massive enough to pro- 
duce these shells would have been easily detected in 
our data. We suggest that the supershells are indeed 
formed internally, but that some additional energy- 
boost (e.g. through magnetic fields) is needed to ex- 
plain their large energies. 
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